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ABSTRACT 

We examine the possibility of evolution with redshift in the mean rest-frame 
ultraviolet (UV; A < 4500A) spectrum of Type la Supernovae (SNe la) sampling 
the redshift range 0<^<1.3. We find new evidence for a decrease with redshift 
in the strength of intermediate-mass element (IME) features, particularly Si II 
and to a lesser extent Ca II 'H&K' and Mg II blends, indicating lower IME 
abundances in the higher redshift SNe. A larger fraction of luminous, wider 
light-curve width (higher "stretch" ) SNe la are expected at higher redshift than 
locally, so we compare our observed s pectral evolution w ith that predicted by a 



redshift-evolving stretch distribution (IHowell et al.ll2007l ) coupled with a stretch- 



dependent SN la spectrum. We show that the sense of the spectral evolution 
can be reproduced by this simple model, though the highest redshift events seem 
additionally deficient in Si and Ca. We also examine the mean SN la UV-optical 
colors as a function of redshift, thought to be sensitive to variations in progenitor 
composition. We find that the expected stretch variations are sufficient to explain 
the differences, although improved data at z~0 will enable more precise tests. 
Thus, to the extent possible with the available datasets, our results support the 
continued use of SNe la as standardized candles. 
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Subject headings: supernovae: general — cosmological parameters — ultraviolet: 
general 



Introduction 



Type la Supernovae (SNe la) have emerged as the most practical and immediate path for 
improving our understanding of dark energy. Techniques have been perfected for locating and 
studying hundreds of events over the redshift range 0<z<1.5 where dark energy manifests 
itself, and ambitious new facilities are being planned that will take this progress to the next 
level of precision. However, it remains the case that the use of SNe la as cosmological probes 
is purely empirical, exploiting relationships between photometric properties such as peak 
l uminosity, color and light curve width, that result in distance estimates precise to ~ 7% 
(fairy et ali 120071 : l.lha et all [20071 : IConlev et al1l2008h . 



Early work by lHamuy et al.l (119951 . Il996l ) demonstrated that some of these photometric 
properties depend on their host galaxy type, with SNe la in spirals possessing wider light 
curves (higher "stretch") and are more luminous than those in older elliptical hosts. This 
has subsequently been categorized in te r ms of (perhaps distinct) populations of 'prompt' and 
'delayed' SNe la (IMannucci et al.ll2005l . 120061 : IScannapieco fc Bildstenl 120051 1 . whose proper- 
ties appear to be governed by the sta r formation activity of the host ga laxy and therefore age 
of the progenitor stellar population (jHoweilll200ll : ISullivan et al.ll2006l ). With star formation 
increasing at higher redshift, these dem ographic variations p roduce a shift to mo re luminous, 
wider light curve events with redshift (jHowell et al.l 120071 ). ISarkar et al.l (120081 ) discuss the 
implications and possible biases introduced in dark energy studies. 

A less well-understood variation is the scatter in the rest-frame SN la ultra yiolet (UV) 
spectra from one event t o another. Althou gh model predictions differ in detail (e.g. lHoflich et al 
19981 : iLentz et al. l hoool : ISauer et al. l boosl ) , the form of the UV spectrum short-ward of 4000A 
is believed to be sensitive to progenitor composition and explosion physics. Given the dif- 
ficulty of studying this wavelength range in nearby events, the most comp rehensive studies 
have been accomp li shed using UV da ta redshifted into the visible at z ~0.5 (jEllis et al.ll2008l : 



Foley et al.ll2008al ). lEUis et al.l (120081 . hereafter E08) found significant scatter in the UV con- 
tinuum below ~4000A, at a level which exceeds that predicted assuming quite significant 
metallicity variations in the current models. 

As both SN la progenitor age and metallicity likely evolve with redshift, it is natural 
to ask whether any spectral variations likewise appear redshift-dependent. Recent work by 
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Foley et al.l (l2008al . hereafter F08) has tentatively identified a decrease with redshift, at ~2-cr 



significance, in the strength of blend of lines (Fe II, Fe III, Si II) near 4800A in maximum- 
light SN la spectra, which they attribute to Fe III 5129A. In this letter, we examine the rest- 
frame A<4500A spectrum of SNe la over 0<z<1.3 by comb ining data gather ed using Keck at 
z~0.5 (Eo8l) and Hubb l e Spa ce Telescope (HST) at z>\ faiess et alJl2007L hereafter R07). 
Following iHowell et al.l (120071 ). our goal is to investigate whether any spectral variations are 
present that represent an evolutionary trend, and whether these can be explained by the 
observed shift to more luminous events at higher redshifts. 



2. SN la Spectral Samples 



We compile spectra from three sources correspondin g to high, intermediate and low 
redshift (z). The intermediate- z spectra are described in |E08| together with all aspects of 
their reduction. T argets were drawn exclusively from the Supernova Legacy Survey (SNLS; 
Astier et al.l 120061 ) . forming a homogeneous sample representative of the SNLS parent sam- 
ple. Host galaxy contamination is removed using photometry of the hosts, with the flux 
calibration corrected using contemporaneous imaging of the SN light curves. This galaxy 
subtraction is a critical step when examining ground-based spectra of high-redshift SNe, and 
an incorrect continuum subtraction can systematically affect line feature strengths and bias 
conclusions drawn from them; the wide wavelength range of SNLS data allow an accurate 
host removal on individual spectra. We degrade the spectra to a resolution of 20A in the 
rest-frame to match that of the high- z HST spectra discussed below (the typical width of 
SN la spectral features is > lOOA, or > 5 resolution elements). 



The z>0.9 spectra are taken from the HST ACS grism campaigns of iRiess et al.l (120041 ) 
and IR071 Host galaxy subtraction is less problematic as the SNe were resolved from their 
host due to the superior resolution of HST. Any host that may lie underneath the SN is 
subtracted by defining narrow "sky" regions above and below the SN trace and interpolating 
to derive the contamination. The resolution of the grism (R = A/AA ~ 200 at 8000A) 
equates to ~ 20A in the rest-frame at 4000A for a z—1 SN, improving for higher- z events. 
We rebin all spectra to the 20A wavelength scale. 



At low- z, samples are heterogeneous in nature. Both the |R07| and lE08l samples were 
found with "rolling searches" and can be considered complete with regard to host galaxy and 
environment; the same is not true for loca l sear ches which fr equently involve t arget ing known 
galaxies. We replace the sample used in |E08| with that of iMatheson et al.l (120081 . hereafter 
M08), which has the benefit of a homogeneous reduction technique, resulting in spectra that 
are close to spectrophotomteric (M08). For each SN we select the closest spectrum to a 
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phase of +3d (the mean of the higher- z samples). Due to the large apparent size of the host 
galaxies, the surface brightness at the SN position is small. The M08 spectra are optical only, 
so at A 5,3600 A, we s upplement the sample with th e space-based UV s pectra of SN1992A 
fcirshner et al.lll993h . SN1981B (branch et alJll983h and SN2001ba (e.g. £ et"aDl2008bh . 
the only three spectroscopically normal SNe la in the appropriate phase range. 



The light curves of all the SNe la are fit using the SiFTO light curve fitter (jConley et al. 



20081 ) to place their photometric parameters on the same system. The distributions in light 
curve phase (r e g; measured in days relative to maximum light in the rest-frame B band 
divided by light curve stretch), z, stretch (s) and rest-frame B — V color (C), are shown 
in Figure [H The high-z HST spectra are typically at a later phase due to the scheduling 
and search strategy used. Note that the SNLS, a rolling search, produces more early SN la 
spectra than more traditional search methods used at low-z. We restrict our analysis to a 
comparison sample of spectra with -l<r c ff<8 (important given the evolution in SN spectra 
with light curve phase), a stretch range of 0.8<s<1.30, and a color range of -0.2<C<0.35. 
We onl y consider "spectroscop ically normal" SNe, excluding obvious, low-Si SN1991T-like 
events ( jFilippenko et al.lll992l ). The mean properties are given in Table [U 



The stretch distributions are marginally different, with evidence for an inc rease in mean 
stretch with redshift at ~ 1.5cr significance, as expected following the analysis of iHowell et al. 



( 120071 ) . The C distributions of all three samples are consistent, arguing against the stretch 
differences arising primarily from selection effects related to search biases, as these would 
also be manifest in the color histograms. The sample of \ow-z SNe with UV spectra has a 
particularly low mean stretch, with two of these SNe at the lower end of the stretch range 
(SN1992a at 0.82 and SN1981b at 0.89). 



Table 1. Mean properties of the SN la samples 



z 


N 


T cS 


s 


C 


1.16 (0.97 to 1.39) 


12 


3.5 


1.087(23) 


0.007(18) 


0.48 (0.12 to 0.73) 


18 


2.0 


1.042(20) 


0.000(20) 


0.02 (0.00 to 0.04) 


18 


2.8 


0.993(25) 


0.042(28) 


0.01 (UV spectra only) 


3 


4.1 


0.91(06) 


-0.01(0.05) 
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z>0.9 (R07) 
z<0.05 (M08+UV) 
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Fig. 1. — The distribution of various properties of the spectral sample. Upper left: The 
phase distribution of the entire sample. The vertical lines show the limits used to define the 
comparison samples. Upper right: The redshift distribution of the entire sample. Lower left: 
The SN stretch distribution for just the comparison sample. The dotted/dashed/dot-dash 
vertical lines show the weighted means for the SNLS/low-^/HST samples. Lower right: The 
SN color distributions, weighted mean colors marked. 
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3. Spectral Comparisons 



We construct mean spectra for the three samples following lE08l . All spect ra are cor- 



recte d for Milky Way extinction and color-corrected using the SALT-2 color law (IGuy et al. 



20071 ). reduci ng, th ough not eliminating, dispersion among SN la spectra and aligning them 
in the mean (IE08I ) . The spectra are normalized to have the same flux through a box fil- 
ter defined from rest-frame 3750 to 4100A (denoted IP), chosen as a common wavelength 
region for all spectra (changing this wavelength range does not change our results). The 
error on the mean spec trum is estimated using a bootstrap resampling technique (see IE08 . 
for more details). IR07I previously presented a m ean spectrum for the high- z HST spectra 
(used in comparisons with lower-redshift by lF08l ). calculated over a wider phase range with 
different error estimat ion t echniques, and lacking color-corrections on the indiv idual spectra. 
A comparison of this R07 mean with (e.g.) the intermediate-^ mean of E08l would not be 
appropriate, so we rederive the HST mean for comparison on a consistent basig^- 

The comparison of the mean spectra is shown in Figure [2J As noted by previous au- 
thors conducting similar analyses, the spectra appear very similar. In particular, with the 
exception of the \ow-z UV mean spectrum, the agreement in the shape of the continua is 
remarkable, suggesting the color corrections made to the spectra based on their optical color 
are effective in aligning the spectra across the entire wavelength and z range. However, we 
find new evidence for variations in some ion feature strengths, most obviously around 4000A 
at the position of the blueshifted Si II 4128A feature. Since the luminosity of the SN la is 
controlled by the amount of 56 Ni synthesized in the explosion, for a Chandrasekhar-mass ex- 
plosion, brighter explosions imply increased 56 Ni production at the expense of intermediate 
mass elements (IMEs)E such as Si, Ca and Mg. So to first order one would expect to see 
fewer IMEs in more luminous, or broader- light curve SNe. Indeed, the Si II 4128 feature has 
been shown to h ave a lower equivalent width in broader-lightcurve SNe la (IGuv et al.l 12 007; 



Ellis et al.ll2008l ) and, by extension, in SNe la originating in spiral galaxies (IBronder et al. 



2008). Given the observed drift in mean stretch with z of the SN la population (IHowell et al. 



20071 ). arising from shifting demographics of the SN la population, a generic prediction of 
two-component prompt+delayed SN la rate models is that the amount of IMEs in genera l, 
and Si in particular, in the average high- z spectrum should be smaller (IHowell et al.l 120071 ) . 
This is in the same sense as the trends observed in the comparison of the means in Figure [2j 



lr The |EQ3 SNe used are: HST05Str, HST05Gab, HST05Lan, HST04Pat, HST04Gre, HST04Eag, 
HST04Sas, SN2003az(Tor), SN2002fw(Aph), SN2003dy(Bor), HST04Omb, sn2002dd. 

2 Stable, nuclear statistica l equilibrium eleme nts like 54 Fe or 58 Ni play a role, but they are relatively 
constant from SN to SN (e.g. iMazzali et alj|20071 ). and have only a small effect on this argument. 
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Fig. 2. — The comparison of the Keck-SNLS (black), low-z (blue) and HST (red) mean 
SN la spectra. The resolution is 20A and no additional smoo thing is performe d. The error 
bars are derived using a bootstrap resampling technique (see Ellis et al.l I2008L for details). 
At A < 3500A in the low-z SNe, no error bars are plotted as only three SNe enter the sample. 
The vertical arrows show the boundaries of the V (right) and £702' (left) filters. 
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The interpretation of SN la spectra is complicated, and this simple picture may not tell 
the whole story. Observing a spectral feature in absorption only constrains the amount of 
material in a given ionization state, along the line of sight, above the photosphere at the time 
of the spectrum. At maximum light , approximately 2/ 3 of the mass of the SN la remains 
hidden under the photosphere (e.g. iHoflich et al.lll998l ). One can imagine a pathological 
situation in which more luminous SNe la have fewer IMEs above the photosphere, and 
more below, but this would be at odds with the basic physical principle that fusion in denser 
regions produces heavier elements, a generic property of SN la models. Instead, a more likely 
complication is that more luminous SNe are hotter, and thus more ionized. In this situation, 
in the hottest regions Si II may be ionized to Si III, resulting in lower Si II equivalent widths 
(EWs). Fortunately, this potential complication does not change the basic principle that 
lower EWs of singly ionized IMEs are expected in the spectra of more luminous SNe la. 



F08l . comparing z=0 and ,2=0.2 -0.8 mean spectra, find a lower strength in the Fe II 
blend at 4800A in higher- z spectra. IF08I state they are unable to distinguish evolving SN 
la demographics with evolution in the SN properties. At first sight, their trend appears 
consistent with the results of this paper: higher-stretch SNe la will be hotter, and thus Fe 
ionized to Fe III at the expense of Fe II. However, the 4800A feature is a comple x one, a blend 



of Si II (5041A J3055A), Fe II (4923A, 5018A, 5169A) and Fe III (5129A) (e.g. lHatano et al. 



19991 ). In fact, lF08l attribute the weakening strength of the red side of the feature to Fe III 
which would be inconsistent with the simple picture outlined above. Given the presence 
of the other Fe II and Si II lines in this 4800A feature, further detailed modeling work is 
required before a robust conclusion with regards to this feature can be reached. 

We investigate the variations in more detail by measuring the strength of three accessible 
IME spectral features (the blend of Ca II at rest-frame ~3950A and Si II at ~3860A, Si II at 
4128A, and Mg II at 448lA). We fit a Gaussian on a lin ear pseudo-continuum background and 
calculate the (pseudo) EW (e.g. iGaravini et al.ll2007l ). The error is estimated by performing 
the EW fit on each of the bootstrap resampled spectra and taking the standard deviation of 
the resulting distribution. The variation with redshift of these features is shown in Figure (3) 
the trends are for a decreasing EW for all three ions with increasing redshift. We emphasize 
that these measurements are made on mean spectra; in many cases the individual high-z 
spectra show evidence of Si II (1R071 ). 



Further into the UV, the continuum level shows a redshift-dependent trend, with the 
spectra broadly consistent from 0.5<z<l, but depressed at z<0.05. Potentially this is an 
important result; however the different SN la samples possess different mean stretches, and 
the UV-optical color of the high er-stretch samples may be expected to appear bluer (e.g. 
Guv et al.ll200~7l : lEllis et alJlioOa l. Furthermore the paucity of local data and the intrinsic 
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Fig. 3. — The variation with redshift of various features in the three mean spectra. Top 
panel shows the L r 02'-L r ' color measured from the spectra, and the lower three panels show 
the equivalent width (EW) of the three features Ca II, Si II and Mg II. Redshift errors are the 
standard deviations of each sample. Overplotted in red is the predicted variation from the 
observed shift in SN demographics (stretch) with redshift, normalized to the intermediate-^ 
sample; the hashed area shows the uncertainty in this prediction based on the errors in the 
intermediate-^ measurement. 
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UV scatter from one event to another ( E08I ) may explain such discrepancies. We quantify 
this "evolution" by measuring a U V-optical color on the mean spectra with two synthetic 
filters denoted {702 ' (following IConley et al.ll2008l ) and {7', both marked in Figure [2J Errors 
are again estimated by analyzing the bootstrapped spectra, and the results are plotted in 
Figure [31 



4. Discussion 



A key question is whether the changes seen in Figures [2] and [3] are a signature of evo- 
lution in individual SN la properties, or instead a product of changing demographics in the 
SN population, i.e. the relative increase in SNe la from younger progenitors at higher red- 
shift. Are the spectral properties of SNe la, at a fixed stretch, evolving with redshift, or 
is the previously observed redshift drift in the mean stretch causing the observed spectral 
evolution? To address this, we predict how the various spectral features might be expected 
to ev olve with redshift w ere a demographic change solely responsible. We use the analy- 
sis of iHowell et all (120071 ). which predicts how the fraction of the two SN la components 
(e.g., 'A+B', delayed+prompt, or old+young) , and hence the mean SN la stretch, changes 
with redshift. This work is based on the galaxy-type dependent stretch distributions of 
Sullivan et all (120061 ) coupled with a cosmic star formation h i story and various sim ple two- 
component SN la rate models (jScannapieco fc Bildstenll2005l ; iMannucci et al.l 120061 ). 



We construct a mean spectrum for low and high-stretch SNe la based on the sample of 

E08l . defining low -stretch (s < s n ) an d high-stretch (s > s ) sub groups. We use a split point 



of s = 1.0 (e.g. IHowell et al.l 120071 ). giving si ow = 0.92 and Shigh = 1-09. The low-stretch 



mean is assigned to the "A" (delayed) component, and the high-stretch mean to the "B" 
(prompt) component. A model mean spectrum for any redshift can then be constructed by 
combining the two stretch-dependent means in the appropriate ratio of the two components. 
Predictions for the behavior of the various spectral features can be made by measuring 
these features on mean spectra constructed at different redshifts, and are over-plotted on 
Figure [3] (there is some circularity here, as the sample of spectra are also used to make the 
intermediate- z measurement; however here we are interested in the predicted change with 
redshift rather than the absolute values). The sensitivity of the predictions to the stretch 
split point is explored varying sq by 5%; this has a very minor effect on the results. 

The predicted tracks of the model generally match the observations, although the un- 
certainties remain large. The decrease in Si II EW with redshift is reproduced, though there 
is evidence that the high-z mean shows an additional deficit in both Si and Ca compared to 
the lower-redshift spectra. This could also reflect an increase in the fraction of SN1991T-like 
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events, believed to arise from a younger population, in the z>l sample, where excluding 
these objects is more difficult due to the lower signal-to-noise of the spectra. Th e predicted 



change in UV-optical color fails to match the data, though the robust datasets (IE08I ; IR07I ) 
show only modest changes consistent with that implied by the change in mean stretch. The 
main discrepancy occurs with the local data. Such a discontinuous behavior seems unlikely 
and probably arises from the incomplete sample of local UV spectra, re-emphasizing the 
need for local, spectroscopically-normal, SNe la observed in the UV. 

We have shown, for the first time, that the strength of some SN la intermediate-mass 
element spectral features evolve systematically with redshift. However, we hypothesize that 
this evolution is related to the p reviously identified shift in SN la photometric properties 
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with redshift flHowell et all 120071 ). and we show it can be partially model ed using a simple 



two-component SN la rate model and the stretch dependent spectra of IE08I . After accounting 
for these expected trends, some small discrepancies in the highest redshift spectra remain, 
though none at high significance. Thus, we conclude that the observed evolution in SN 
la spectral features can be explained by, and is consistent with, changes in SN la demo- 
graphics, rather than evolution in the properties of individual, photometrically-similar SNe. 
Controlling for these demographic shifts is the challenge for SN la precision cosmology ex- 
periments, for example by ensuring that an appropriate fc-correction template is used as a 
function of redshift, and confirming that the empirical relationships used to calibrate SNe 
la in cosmological applications are valid across all environments. 
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